Fractional Boltzmann equation for 
resonance radiation transport in plasma 
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Abstract: The fractional Boltzmann equation for resonance radiation transport in plasma is 
proposed. We start from the standard Boltzmann equation, averaging over frequencies leads 
to appearance of fractional derivative. This fact is in accordance with the conception of latent 
variables leading to hereditary and non-local dynamics (in particular, fractional dynamics). The 
presence of the fractional material derivative in the equation is concordant with heavy tailed 
distribution of photon path lengths and with spatiotemporal coupling peculiar to the process. 
We discuss some methods of solution of the obtained equation and demonstrate numerical results 
in some simple cases. 
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1. INTRODUCTION 

The standard (Markovian) transport model based on 
the Boltzmann equation can not describe some non- 
equilibrium processes called anomalous that take place 
in turbulent plasma, interstellar magnetic fields, disor- 
dered semiconductors, and other complex structures (see 
Bouchaud k Georges (1990); Metzler & Klafter (2000)). 
Causes of anomality lie in non-uniformly scaled (fractal) 
spatial heterogeneities, in which particle trajectories take 
cluster form. Furthermore, particles can be located in some 
domains of small sizes (traps) for a long time. Estimations 
show that path length and waiting time distributions are 
often characterized by heavy tails of the power law type. 
This behavior allows to introduce time and space deriva- 
tives of fractional orders. Distinction of path length distri- 
bution from exponential is interpreted as a consequence of 
media fractality, and analogous property of waiting time 
distribution as a presence of memory. 

During last decades, essential progress in description of 
the asymptotical behavior of such processes is achieved by 
passage from the standard diffusion equations to the equa- 
tions with fractional derivatives (see, for example, Met- 
zler & Klafter (2000), Samko et al. (1993), Hilfer (2000), 
Podlubny (1999), West et al. (2002), Nigmatullin (2006), 
Sibatov & Uchaikin (2009), Uchaikin et al. (2008), Datsko 
& Gafiychuk (2010)). The diffusion equation describes 
only asymptotical part of the solution of the Boltzmann 
equation. For more complete description the passage to 
fractional derivatives should be performed not from the 
diffusion equation but on the previous stage, i. e. from the 
Boltzmann equation. 

Nonnenmacher & Nonnenmacher (1989) proposed frac- 
tional generalization of the bilinear Boltzmann equation 
in the form 
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where /(r, v, £) is the one-particle distribution function, 
C(/, /) is the collision integral, and rjD t " is the Riemann- 
Liouville partial derivative of fractional order 1 — v, < 
v < 1. The derivation of this equation stated in Non- 
nenmacher & Nonnenmacher (1989) is rather formal. It is 
based on the analogy between corresponding standard dif- 
fusion equation and its fractional generalization. Neverthe- 
less, it leads to valid asymptotic equations whose solutions 
are consistent with Monte Carlo simulation. However, one 
should clearly understand what the equation describes. 
Authors Nonnenmacher & Nonnenmacher (1989) have not 
paid sufficient attention to this issue. The fractional equa- 
tion describes hereditary kinetics with memory kernel of 
power law type which is usually caused by waiting times in 
local domains with power law asymptotes in distribution 
function. The process does not assume spatial-temporal 
coupling: the waiting times are considered to be indepen- 
dent of random free path lengths distributed according 
to exponential law. Moreover, the detailed consideration 
concerned only one-dimensional process. 

The linear Boltzmann equation describing the anomalous 
transport of neutrons within processes of scattering, fis- 
sion, and absorption, was generalized to fractional form by 
Kadem & Baleanu (2010) and its solution was considered 
under more realistic conditions. 

We are going to consider fractional generalized kinetics 
in connection with radiation transport in plasmas. The 
standard Boltzmann equation has found very important 
applications in studying transport of particles in plasma. 
It is known (see Pereira et al. (2004), Berberan-Santos et 
al. (2006), Pereira et al. (2007)) that a stochastic behavior 
of atomic resonance radiation migration in plasma is of 



superdiffusive type. Pereira ct al. (2004) have shown that 
photon trajectories for Doppler, Lorentz, and Voight line 
shapes under the assumption about complete frequency 
redistribution are Levy flights. Superdiffusivity is a conse- 
quence of the heavy-tailed nature of path length distribu- 
tion 

P(0 = Prob{i? > £} oc f. 

The lighter the tail of the line shape function, the less 
heavy the tail of path length pdf. Berberan-Santos et al. 
(2006) concluded that this tail is always heavy [v < 1). 

Finiteness of photon velocities assumes that the process 
is Levy walks, not Levy flights. Levy walks give a proper 
dynamical description of the superdiffusive behavior. The 
temporal and spatial variables of Levy walks are strongly 
correlated. Zolotarev et al. (1999), Zaburdaev & Chukbar 

(2002) , Sokolov & Metzler (2003), Chukbar & Zaburdaev 

(2003) , Uchaikin & Sibatov (2004), Uchaikin & Sibatov 
(2009)) proposed equations for Levy walks, but these 
equations are obtained in diffusion approximation and 
authors concentrate on the one-dimensional case. 

In the present paper, we derive the fractional Boltzmann 
equation for propagation of resonance radiation in plasma. 
We start from the standard Boltzmann equation, which 
suppose exponential distribution of path lengths of photon 
with a given frequency uj. Averaging over the random 
frequency leads to power law distribution of path length 
(see Pereira et al. (2004)) and to appearance of fractional 
derivative in the equation. This fact is in accordance with 
the conception of latent variables leading to fractional 
dynamics (see Uchaikin (2008)). It is interesting that equa- 
tion contains three-dimensional material derivative gener- 
alizing the one-dimensional case introduced in Sokolov & 
Metzler (2003). Finally, we discuss methods of solution of 
the obtained equation and present results in some simple 
cases. 



2. PECULIARITIES OF THE PROCESS 

The process under consideration looks as follows. Initial 
photon with frequency uj' and velocity c entering the gas 
goes a random free path R distributed according to the 
exponential law 

Prob{i? > = exp(-o-(w)£)- 
Then, it interacts with an atom which absorbs the photon 
passing to excitement state. After a random time T, 
the atom re-radiates a photon with a frequency to being 
close to the initial uj but not coinciding with them. The 
irradiated photon continues its way as before. 

Therefore, the first special feature of the process is a finite 
speed of light. The second one is the random delay of 
the photon rc-cmission. And the third, most important 
feature of resonance excitation transport in a spectrum 
line is the fact that frequency dependence of absorption 
factor <t(oj) repeats a shape of the very radiation spectral 
line w(oj). This leads to a specific situation: most of 
irradiated photons have very short free paths and they 
are absorbed again, practically at once, so that only small 
fraction of photons having frequencies in tails of spectrum 
(far from the center of the line) can go far from re-emission 
points. We observe Levy walk. 



The resonance radiation transport theory develops in two 
directions - in light-engineering and in astrophysics. In the 
first of them, the finiteness of speed of light is neglected, 
that is, one accepts c = oo. The main equation in this case 
is the Biberman-Holstcin equation. The second direction 
concerns astrophysical scales, so that the traveling times 
are much larger then the excitation times and the latter 
can be neglected. Here we concentrate our attention on 
this last problem. 

3. INTEGRAL TRANSPORT EQUATION 

For the sake of convenience, they usually identify radiated 
and absorbed photons, in other words, they say that an 
atom absorbing a photon radiates " the same photon" with 
frequency and direction changed after an interaction with 
the atom, or shorter, the photon is scattered on the atom. 
This process is described by the linear Boltzmann equation 
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which can be transformed to a pure integral form: 
$(r,n,w,i) = $o(r,n,w,i) 
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xa(uj')<S>(r - SI £, SI', J, t - i/c). 
Here $(r, SI, uj, t)drdSlduj is the mean number of photons 
at the moment t in volume dr — dxdydz flying into solid 
angle dSl with frequencies in [uj,uj + duj) Further, we use 
the standard assumptions. 

• Assumption 1. The complete frequency redistribution 
takes place: 

w(uj <— uj) — w(uj) 

• Assumption 2. The local thermodynamic equilibrium 
is valid: 

cr(uj) = Cw(uj) 

Passing to the photon collision density 

F(r, SI, uj, t) = ca(uj)$(r, SI, uj, t), 

we arrive at 

F(r,Sl,uj,t) = F (r,Sl,uj,t)+ 

OO 

+ J d£ w(uj)cr(uj)e-' T ^ J dSl' J duj' W(Sl <- Sl')x 
o 

xF(r- SI £,Sl',oj',t-t/c) 
Integrating over frequency 

oo 

F(r,Sl,t) = J F(r,Sl,uj,t)duj, 
o 

we obtain the following integral equation 
F(r,n,t)-F (r,n,i) = 

OO 

= J d£ p(0 J dSl'W(Sl <- n') F(r -Sl£,Sl',t- £/c), 



where the kernel is given by 



p(0 = (a(Lu)e- a ^^ = J w(u;)a(uj)e- a ^duj. 



By integrating by parts, the last operator can be reduced 
to the form 



Characteristic feature of resonance radiation propagation 
is the fact that absorption spectrum (as an approximation 
of thermodynamic equilibrium and total redistribution of 
frequency) is proportional to the spectrum of radiation 
a(w) = Cw{uj) 

Calculations lead us to the following asymptotic properties 
of the path length distribution 
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To satisfy this condition, we choose the distribution in the 
form of " fractional exponents" : 

P(0 = 7£"- 1 ^(-70, (3) 

where E a ^(x) is the two-parameter Mittag-LefHer func- 
tion (see, for example Uchaikin et al. (2008)). 

4. FRACTIONAL BOLTZMANN EQUATION 

Generalized Boltzmann equation for the case of path 
length distribution (3) is obtained in the form 

(d \ v 
-+ C nVj $(r,n,i)+/i$(r,n,t) = 



$o(r,n) S(t) 
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where fi = cy and the operator 
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is the fractional generalization of the material derivative. 

Ordinary material derivative is determined as 

f(r,t)-f(r-ch,t-h) 
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The operator 

d_ 

dt dr y 

is the infinitesimal operator generating the semigroup of 
operators 

T h f(r,t) = f(r-ch,t-h) 
Consequently, we have for (— A) a 

(-A) a f(r,t) = (l+c^-Y f(r,t) = 



ds 
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If we consider the case of a space-homogeneous distribu- 
tion of particles $(r, fi, t) = $(fi, t) this operator becomes 
the fractional Rieman-Liouville derivative of order v 

Q + cnv) $(r,n,t)H> D t v $(n,t). 

In a stationary problem, when the function $(r, fi) docs 
not depend on time, 



+ cOV) $(r,fl)4 c(fiV)^(r,fi,(), 



this operator represents fractional generalization of the 
directional derivative. 

5. SOLUTIONS OF FRACTIONAL BOLTZMANN 
EQUATION 

5.1 Analytical results 

In one-dimensional case, when photons can fly only along 
of x-axis, the fractional generalization of the Boltzmann 
equation is of the form 

d d \ v 

_ + c o x _j <f>(x,{i x ,t) + fj,$(x,n x ,t) = 



d d 
= fiJ(x,il x ,t) + [ — + cQ x — 



S(x, £l x , t). 



Using expressions for the collision integral 

J(x,n x ,t) = i[$+(M) + *_(a;,t)] [5(n x -l)+S(Q x + l)} 

and for the onc-dimcnsionally "isotropic" instantaneous 
source 

S(x, Sl x , t) = J [S(Q X - 1) + 5{Q X + l)]S(x)d(t) 



we obtain the following equations: 
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From last equations one follows 
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In the asymptotics of large times the last equation passes 
into the following one 



d 
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solution of which is expressed in terms of elementary 
functions (see Uchaikin & Sibatov (2009)) 

, . , 2 sin Tiv 
= x 

(l-x 2 /c 2 t 2 Y^ 
X (1 - xjct) 2v + (1 + xjct) 2v + 2(l-x 2 /c 2 t 2 ) v COS TTV ' 

Quantitative analysis of distributions can be performed 
using the method of moments. Due to finiteness of photon 
velocity, its position is bounded random variable and 
according to Cramer's criterion distribution of excitations 
is unambiguously determined by the moments of this 
distribution. Moments can be found from a characteristic 
function. 

For the one-dimensional case, the Fourier-Laplace trans- 
formation of the equation (6) leads to the following ex- 
pression 

(A - ickf- 1 + (A + icky- 1 + — (A 2 + c 2 k 2 Y 

A) = ^ . 

(A - icky + (A + icky + — (A 2 + c 2 k 2 Y 

A* 

In the asymptotics of small k and A corresponding to the 
diffusion limit, we obtain 

~ (A-icfc)"- 1 + (A + i c fc)"- 1 

1 ' ' ~ (a - icky + (a + icky ' 

Moments are determined by the relation 

_d"d>(fc,A) 



m n (A) 



dk n 



In isotropic case, Fourier-Laplace transformation of (4) 
leads to the following asymptotical form for <&(k, t) 

J (A - icfikf -1 dil 



*(M) 



5.2 Numerical results 



J (X-icSlkfdM 



Here, a computational algorithm for the process is de- 
scribed, realized in the Monte Carlo code and used for 
computing the resonance radiation and excitation distri- 
bution in plasma. An important advantage of the method 
is simplicity of taking into account boundary conditions 
and inhomogeneity of the medium. 

This method represents randomized scheme of realization 
of generation expansion. It can be considered as statistical 
simulation of transport process characterized by known 
distributions of independent random parameters determin- 
ing a random trajectory. Set of large number of indepen- 
dent trajectories allows to estimate statistically a process 
quality interesting for us. 




x/vt 



Fig. 1. Comparison of analytical solutions with the re- 
sults of Monte Carlo simulation for isotropically one- 
dimensional case as functions of dimensionless vari- 
able x/ct. 
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Fig. 2. Results of Monte Carlo simulation for isotropic 3D 
case. 

There is the typical sequence of operations for simulation 
of photon trajectory. 

— A position of excited atom at the initial time moment 
is chosen. PDF of its radius- vector Ro is determined by 
given initial distribution of excitations iV(r, 0) through 
the relation p R = N(r, 0)/ / N(r, 0)dr. 

— The random time till the moment of photon emission 
by the atom is generated according to the exponential 
distribution pr(t) — 7e _7t . 

— The frequency of emitted photon is distributed with the 
pdf w(oj). 

— A direction H of photon emission is generated. This 
direction is characterized by random variables cos and 

which are uniformly distributed (in isotropic case) in 
[—1,1] and [0,27r), respectively. 

— A photon path length R is simulated in infinite homoge- 
neous medium, pr(£) = a{oj)e~ a ^ £ - . If the whole seg- 
ment [Rq,Ro + OiJ] keeps within the volume occupied 



by plasma, its second end is taken as the position of 
another excited atom. This atom absorbs the photon. If 
this segment traverses absorbing or transparent bound- 
ary, the trajectory comes to the end. 
— For new excited atom, waiting time till the moment 
of photon emission is simulated. And so on, described 
operations are repeated before the end of trajectory. 

From physical point of view, the fractional operator in 
the Boltzmann equation means that anomalously long free 
paths arising from time to time tear the trajectory into 
more or less localised clusters which look separated at any 
scale. 

Solutions obtained by Monte Carlo simulations are pre- 
sented in Fig. 1, 2 (points). One-dimensional solutions 
listed in the previous subsection agree with the numerical 
results. 

6. CONCLUSION 

The fractional generalization of the Boltzmann equation 
taking into account spatiotemporal coupling peculiar to 
resonance radiation transport is proposed. This coupling 
is provided by heavy tailed distribution of photon path 
lengths which arises after averaging over frequencies and 
by finiteness of photon velocity. Heavy tails of path length 
distribution leads to appearance of fractional derivative in 
the kinetic equation, but due to coupling, this fractional 
operator represents not simple fractional Laplacian. It is a 
fractional analogue of the material derivative that is agree 
with results obtained by Sokolov & Metzler (2003) for one- 
dimensional Levy walks. 
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